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基于 小 波 变换 的 卫星 阻力 系数 分 析 


刘 P, XI, WAA, MS, EE 
( 中 国 科学 院 空间 科学 与 应 用 研究 中 心 ， 北 京 100190) 


摘 要: 以 不 同 弧 长 解 算 天 宫 一 号 空间 实验 室 连 续 55 天 ( d) GPS 轨道 数据 得 到 卫星 阻力 系数 C, 序列 ,并 对 
同期 的 太阳 辐射 通 量 fi, ;和 地 磁 活 动 指 数 a, ,4, 序列 高 、 低 频 多 层 小 波 分 解 。 空 间 环境 参量 和 C, 序列 各 层 信号 滑 
动 相关 性 的 统计 分 析 结 果 表 明 小 波 分 解 的 阻力 系数 序列 与 对 应 的 fi,;,4, 序列 整体 有 具有 很 好 的 滞后 相关 性 ,验证 
了 通过 各 ,4, 预报 C, 的 可 行 性 。 对 不 同 回归 模型 和 回归 分 析 样 本 序列 长 度 n, 分 析 C, 预报 值 残 差 表明 : 采用 fo， 
一 元 回归 模型 ,并 取 n=20 dif, C, 预报 结果 最 优 。 另 外 , NRLMSIS - 00 模型 在 350 km 高 度 的 中 低 纬度 区 域 ,地 
磁 活动 上 升 相 响应 过 度 ,而 下 降 相 则 响应 不 足 ; 太阳 辐射 通 量 上 升 相 响应 适度 ,下 降 相 响应 过 度 。 最 后 对 应 用 前 景 
进行 总 结 ,该 方法 对 提高 航天 器 轨道 预报 精度 具有 重要 参考 价值 。 
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Analysis of Satellite Drag Coefficient Based on Wavelet Transformation 


LIU Wei, WANG Rongdan,; LIU Si-qing, SHI Li-qin; GONG Jian-eun 
( Center for Space Science and Applied Research, Chinese Academy of Sciences, Beijing 100190, China) 


Abstract: Drag coefficient sequence C, is obtained by solving Tiangong 1 continuous 55 days ( d) GPS orbit data with 
different arc length. The same period solar radiation flux fi, ; and geomagnetic index a,; A, sequences are high and low 
frequency multiowavelet decomposed. Statistical analysis results of the multilayer sliding correlation. between space 
environmental parameters and decomposition of C, show that lag correlation between the satellite drag coefficient sequence 
after wavelet decomposition and the corresponding level of fj, ;, A, sequence is good. lt is also verified that the C, 
prediction is feasible. Prediction residuals of C, with different regression models and different sample lengths are analyzed. 
The results show that the case setting sample length of 20 days and selecting fj; regression model is best. It is also show 
that NRLMSIS-00 model s response in the region of 350 km ( Tiangong' altitude) and low-middle latitude ( Tiangong s 
inclination) is excessive during the ascent phase of geomagnetic activity A, and is inadequate during decline phase. 
Additionally, for the low-frequency decomposition sequence NRLMSIS-00 model$ response is appropriate during the ascent 
rising phase of fio7. For the high frequency decomposition sequence; NRLMSIS-00 model's response is small-scale 
inadequate during the ascent phase of fio; and is excessive during the decline phase of fio7. Finally; the summary of 
potential use and outlook are listed; this method has an important reference value for improving the spacecraft orbit 
prediction accuracy. 
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"ET 的 主要 误差 源 。 目 前 典型 大 气 模 式 中 误差 均 在 
15% ~ 3096  -? ,太阳 活动 峰 年 或 磁 暴 期 会 更 大 。 
大 和 气 阻 力 摄 动 是 近 地 空 间 目标 轨道 确定 、 预 报 该 误差 是 经 验 模式 所 固有 ,由 空间 环境 指数 建 模 不 
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完善 引起 , 可 通过 解 算 同期 卫星 轨道 数据 补偿 修 
正 "”。 大 气 模式 短期 修正 方面 ,国内 外 学 者 做 了 大 
量 研究 工作 。Doornbos 等 利用 TLE 轨道 数据 中 丰 
富 的 阻力 信息 ,对 CIRA - 72 模式 进行 短期 修正 ,并 
使 用 5 颗 近 地 点 280 -530 km 目标 验证 修正 效果 ， 
表明 仅 一 天 的 参数 修正 可 减 小 密度 中 误差 15% ~ 
30% ,考虑 其 它 参 数 的 影响 中 误差 至 少 能 降低 
12% Ë, Lechtenberg 等 基于 精密 星 历 对 短期 大 气 
模型 密度 进行 修正 ”。Pardini 等 通过 第 23 太阳 活 
动 周 峰 年 4 颗 卫 星 半 长 轴 变 化 ,分 析 了 6 种 (JR - 
71, MSISE-90/00 , GOST -2004 和 JB06/08) 大 气 模 
式 偏差 ,指出 现 有 经 典 模式 均 高 估 500 km 以 下 大 气 
密度 约 7% - 2096 P ,此 偏差 是 建 模 过 程 中 考虑 卫 
星 阻力 系数 高 度 无 关 引 起 "3 。 同 时 文献 5 -8] 
指出 卫星 阻力 系数 随 轨道 高 度 、 卫 星 构 形 和 材质 等 
因素 变化 。Pardini 等 分 析 轨 道 高 度 150 ~ 1500 km 
之 间 11 颗 球 形 卫星 整个 太阳 活动 周 轨道 衰减 数据 ， 
评估 JR -71 和 MSISE - 90 模式 , 认为 两 模式 误差 
空间 环境 变化 及 轨道 高 度 是 相关 的 , 且 任何 空间 
环境 下 , 两 模式 均 高 估 400 km 以 下 大 气 密 度 晶 。 
大 气 密 度 除 包含 11 年 .27 天 、 半 周年 、 周 年 半 
周 日 、 周 日 和 潮汐 等 错综复杂 的 周期 变化 外 ,还 存在 
磁 暴 等 因素 引起 的 非 周 期 变化 ,有 些 因素 难以 模 制 ， 
使 得 大 气 密度 精确 建 模 工作 相当 困难 。 文 献 [0 - 
11] 指出 现 有 大 气 模式 误差 与 空间 环境 参量 相关 ， 
并 建立 基于 远 柴 外 能 量 EQ. 指数 、 地 磁 指 数 a, 的 短 
期 动态 密度 修正 系数 。 大 气 阻 尼 摄 动力 计算 式 中 售 
大 气 密 度 、 相 对 速度 和 弹道 系数 等 因素 。 研 究 表明 ， 
一 定 弧 段 内 的 大 气 密度 误差 可 通过 阻力 系数 Cv F 
滑 补偿 ,以 削弱 大 气 阻 力 误差 。 陈 建 荣 等 在 磁 暴 期 
间 用 折线 型 C 吸收 大 气 模式 误差 ,提高 了 轨道 确定 
Km "^U 。 刘 舒 蔡 等 分 析 磁 静 期 阻力 系数 与 地 磁 
指数 和 治 迹 残 差 的 关系 ,建立 中 长 期 轨道 预报 中 阻 
力 系数 的 补偿 算法 "“ 。 
小 波 变 换 源 于 20 世纪 80 年 代 地 球 物理 研究 中 
地 震 信 号 分 析 ,其 克服 傅立叶 变换 的 缺点 ,具有 
时 域 和 频 域 局 部 化 特征 ,变换 后 可 获取 目标 信号 各 
频率 子 段 信 息 。 小 波 变 换 是 强 有 力 的 时 频 分 析 工 
上 县 ,有 人 将 其 称 为 “数学 显微镜 ”。 小 波 分 析 已 广泛 
用 于 信号 、 图 像 处 理 、 模 式 识 别 及 大 气 灌流 等 诸多 领 
域 。 文 献 l6] 利用 小 波 变换 时 域 特性 ,对 降水 时 间 


序列 进行 多 时 间 尺 度 小 波 分 析 , 为 旱 涝 预 报 提 供 了 
重要 信息 。 

天 宫 一 号 ( Tiangong 1, TG1) 是 我 国 首 个 空间 实 
验 室 。 自 2011 年 发 射 以 来 积累 了 大 量 探测 数据 ,为 
空间 科学 研究 贷 定 了 基础 。 基 于 大 气 模式 误差 与 空 
间 环 境 参量 相关 , 且 可 通过 阻力 系数 进行 补偿 的 理 
论 。 本 文 以 不 同 弧 长 解 算 2013 -2 -5 T00: 00 至 
2013 -3 -31 T00: 00( BJT) 间 TG1 连续 55 天 ( d) 
GPS 轨道 数据 ,得 到 卫星 阻力 系数 C, 序列 ,利用 30 
法 则 剔除 异常 数据 。 对 同期 太阳 辐射 通 量 ,名 ， 和 地 
磁 指数 o, ,4 序列 进行 高 .低频 5 层 小 波 分 解 ,分 解 
后 的 各 层 信 号 进行 了 滑动 相关 性 分 析 。 结 果 表明 ， 
小 波 分 解 后 的 C, 与 名 ,4, 各 层次 信号 有 具有 很 好 的 
滞后 相关 性 ,验证 了 C, 通过 /;,;,4, 预报 的 可 行 性 。 
通过 回归 分 析 建 立 C, 与 /io;,4, 关系 ,并 使 用 GPS 
轨道 数据 解 算 C, 值 验证 预报 的 有 效 性 。 最 后 对 方 
法 适用 范围 与 局 限 性 等 进行 讨论 ,该 方法 对 提高 航 
天 器 轨道 预报 精度 具有 重要 意义 。 


1 阻力 系数 解 算 


卫星 在 轨 运 行 中 ,大 气 阻力 加 速度 cu. 可 表达 
为 : 


1 


A drag Tis 5 Bpr; ( 1) 


式 中 : 弹道 系数 B = CAM , CL 为 卫星 阻力 系数 ， 
A/M 是 面积 质量 比 , p 是 卫星 所 处 位 置 大 气 密度 , v, 
是 卫星 相对 大 气 的 速度 矢量 。 

大 气 阻力 加 速度 由 C, A/M , v, 和 p 四 因素 共 
同 构成 。 大 气 阻 力 摄 动 误差 为 各 因素 误差 与 因素 间 
耦合 误差 之 和 。 故 大 气 密度 误差 通过 阻力 系数 补偿 
的 原理 可 表达 为 : 

B[ pria m B pra (2) 
式 中 : 8 -为 轨道 确定 中 估算 弹道 系数 , p 为 大 气 密 
度 模 式 计算 值 , A 为 轨道 拟 合 弧 长 。 

T A/M 已 知 的 情况 下 , 估算 阻力 系数 C, 与 
估算 弹道 系数 B 是 等 效 的 。 地 磁 指 数 a, ,4, 分 为 
3 h 值 和 日 值 ,太阳 辐射 通 量 fi; 为 日 值 , 故 本 文 
分 别 以 3 h 弧 长 和 1d sm A& 82 $3 F8 72 98 7 C, 
序列 。 

1.1 解 算 方法 与 模型 


阻力 系数 解 算 本 质 是 将 其 包含 到 状态 矢量 天 
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中 , 进行 微分 修正 。 通 过 GPS 轨道 数据 中 位 置 矢量 
r, 建 立 测量 方程 : 
= G(i,X) «V (3) 
将 测量 方程 线性 化 ,在 参考 状态 量 X 处 泰勒 
展开 , 略 去 高 阶 项 ,得 到 条 件 方 程 : 
y= Cr+5 (4) 
式 中 : y =r, -六 是 残 差 , r ,r 分 别 为 位 置 矢量 测量 
值 和 计算 值 ; 状态 量 匀 包含 位 置 、 速 度 矢量 和 阻力 
XEGBUX-X(n»C);x-X,-X -X-X 


为 待 估 状 态 量 的 修正 量 ; G 是 相应 m x n 维 高 矩阵 ， 
m 是 个 测量 数据 对 应 的 维 数 , n 则 是 待 估 状 态 量 
的 维 数 ,通常 m 远大 于 n ,充分 利用 测量 数据 统计 
信息 。 
阻力 系数 解 算 中 ,主要 模型 和 参量 : 
1) 70 x 70 阶地 球 引力 场 模型 GGMO2C; 
2) 30 x30 阶 TOPEX 4. 0 海潮 模型 ; 
3) NRLMSIS -00 大 气 模型 ; 
4) 坐标 、 时 间 及 天 文 常数 系统 采用 TERS 
) 
) 


r 


( 2003) 标准 ; 

5) 日 月 引力 摄 动 ; 

6) 太阳 光 压 ,固体 潮 摄 动 。 

通过 多 次 迭代 求解 条 件 方程 ,获取 阻力 系数 最 
优 估 计 , 解 算 连 续 弧 段 得 到 阻力 系数 C, 序列 。 
1.2 异常 数据 剔除 

天 宫 任务 需要 轨道 机 动 或 轨道 衰减 等 , 需 轨道 
机 动 以 保障 任务 或 轨道 维持 。 如 解 算 弧 段 内 含 轨道 
机 动 , 解 算出 的 阻力 系数 就 不 能 准确 表征 大 气 密度 
偏差 , 此 情况 为 阻力 系数 异常 , 表现 为 定 轨 残 差 
RMS( Ru) 较 大 。 通 过 残 差 RMS 统计 ,使 用 3o 
法 则 对 异常 数据 进行 剔除 。 


即 : 如 果 | Reus -aaul| > 3O resina » DUE LER 
该 弧 段 解 算 的 阻力 系数 ,以 保证 获取 的 C, 序列 正确 
反映 大 气 密度 误差 。 其 中 Reina 为 定 轨 残 差 RMS， 
Mesaal 为 Resaa FINI, oua 为 Rasa 序列 标 
E. 
2 数据 处 理 


解 算 阻 力 系 数 分 别 得 到 弧 长 3h 和 1 d 的 Cv 序 
列 , 并 剔除 异常 值 。 对 C, 序列 作 归 一 化 处 理 , 定义 
归 一 化 值 C, = C,/2.2 。 对 C0, 及 同期 /,;,4,,a, 序 


列 进行 高 .低频 5 层 小 波 分 解 ,不 同 尺度 上 分 析 分 解 
信号 相关 性 。 
2.1 离散 小 波 变换 

假设 输入 信号 /(i) , 记 pli) = 27p(274 - 
k) s palt) = 2779(27t - E) > olt) &u(i) 为 多 
分 辨 率 分 析 中 的 尺度 画 数 和 小 波 画 数 ; Ce, (0) 与 
{wi (0) ) 为 正 交 基 夯 数 集 。 记 Pf = /在 第 / 级 上 的 
一 维 离散 小 波 变换 为 ( DWT) 。 通 过 正 交 投影 Pf 与 
QSR P SORR. : 

Pf =P, f+Qf= 2. Cg + 2: di, (5) 


R: d = Yu) dt od = D eln) ih Ut 

PEDEICCIEICON 
为 低 通 和 高 通 滤波 系数 , 由 基 瑟 数 C, (0) 、 
{yalt} 确定 ; p 为 滤波 系数 的 长 度 ; NN 为 输入 信号 
的 长 度 , 工 为 所 需 的 级 数 。 最 终 得 到 各 级 小 波 系 数 
序列 c,d 。 

小 波 基 匡 数 选用 Daubechies, 图 1 为 弧 长 1d 归 
一 化 阻力 系数 C, 及 对 应 4, 序列 多 级 分 解 结果 。 
1( a) 顶部 两 曲线 为 4, 和 C, 原始 序列 , 原 序列 间 不 
存在 显著 相关 性 。 图 1( a,b) 其 余 曲 线 为 两 序列 低 
频 分 解 [5 ~ L1, 低频 部 分 表现 出 一 定 负 相关 性 。 
同样 ,图 1 (c,d) 除 顶 部 曲线 为 A, 和 C, 原始 序列 
外 ,其 余 为 两 序列 高 频 分 解 Hz5 ~ H1, 高 频 部 分 呈现 
出 更 强 负 相关 ,表明 NRLMSIS - 00 模式 在 350 km 
高 度 中 低 纬 度 区 域 ( Tiangong 1, h =350 km, i 
42°) 地 磁 活动 上 升 相 存在 过 响应 ,而 下 降 相 则 响应 
不 足 。 

2 给 出 弧 长 1d BE C, 0 f, 序列 小 波 分 解 
结果 ,图 2( a) 顶部 两 曲线 为 /i。; 和 C, 原始 序列 , 原 
始 序列 间 存 在 较 强 的 相关 性 。 图 2( a,b) 其 它 曲 线 
为 两 序列 低频 分 解 .5 ~ L1. 低频 部 分 L5 太阳 辐射 
AE 上 升 相 C 序列 基本 保持 稳定 , 下降 相 C, 
序列 相应 下 降 , 表明 NRLMSIS - 00 模式 在 350 km 
高 度 中 低 纬度 区 域 , fo, 上 升 相 响应 适度 , 下 降 相 存 
在 过 响应 。 同 样 ,图 2( c,d) H (c) 顶部 两 曲线 为 
f; 和 C, 原始 序列 ,其 它 为 两 序列 高 频 分 解 85 ~ 
Hle foz > 


Daas ;L,k 20513: 


C, 高 频 部 分 呈现 出 很 强 的 正 相 关 , 表明 
NRLMSIS - 00 模式 在 fi, ;上升 相 在 小 尺度 上 的 响应 
不 足 , 而 下 降 相 存在 过 响应 。 
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图 1 地 磁 指数 和 阻力 系数 的 小 波 分 解 图 2 太阳 辐射 指数 和 阻力 系数 的 小 波 分 解 


Fig.1 The multi-wavelet decomposition of A, , C, sequences Fig.2 The multi-wavelet decomposition of fi, ; C, sequences 
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另外 分 析 表 明 , 弧 长 h 解 算 阻力 系数 序列 与 
对 应 w 序列 呈 极 弱 相 关 性 。 可 能 是 解 算 弧 长 不 够 
等 因素 引入 了 其 它 误 差 导致 的 ,限于 篇 幅 具体 分 析 
结果 不 再 列 出 。 

2.2 ”滑动 相关 分 析 

滑动 相关 性 分 析 是 固定 空间 环境 序列 ,前 后 滑动 
C, 序列 并 计算 相关 系数 ,统计 相关 系数 变化 规律 ,为 
C, 预报 分 析 提供 依据 。 首 先 对 4, 及 对 应 的 C, 序列 
作 滑 动 相关 性 分 析 。 选 取 同 一 分 解 层 两 相同 长 度 的 
序列 ,固定 4, 序列 前 后 滑动 C, 序列 , 取 滑 动 值 -3 ~3 
d 计算 两 序列 相关 系数 。 解 算 长 度 55 d 序列 中 , 取 
滑动 序列 长 度 30 4d。 滑动 值 在 -3 ~3 d 范 围 时 ,各 滑 
动 值 均 可 计算 19 个 相关 系数 。 

表 1 列 出 A, , C, 序列 各 层 信 号 不 同 滑动 值 时 ， 
XRB. Rag Rag HS ~ MRa glS ~ Ll 


各 表示 A, , C, 原始 序列 ,高 频 万 ~ 有 ,低频 15 ~ 
L1 序列 相关 系数 均值 。 图 3 是 A, , C, 序列 各 层 分 
解 信 号 滑动 相关 系数 统计 。 图 3( a,g) 是 4,,C, A 
始 序列 滑动 相关 性 统计 , 滑动 值 d 时 原始 序列 间 
表现 出 较 强 的 负 相 关 性 ,如 表 1 所 示 相 关系 数 均值 
-0.3102. 图 3(b ~f) 4418 A, , C, 序列 高 频 ~ 
H1 滑动 相关 系数 统计 ,可 见 原始 序列 及 高 频 #85 ~ 
H2 部 分 滑动 值 取 2 d 时 , 呈现 出 较 强 的 负 相关 性 ， 
相关 系数 均值 分 别 为 : - 0. 3102, - 0. 2322, 
-0. 5576, -0. 37, -0.4396 ,而 高 频 H1 的 相关 性 则 
没有 规律 。 图 3(h -1) 给 出 A, ，C, 序列 低频 [5 ~ 
L1 滑动 相关 系数 统计 情况 ,低频 部 分 情况 类 似 整 体 
表现 出 负 相 关 , 当 滑 动 值 取 2d 时 整体 负 相 关 性 
较 强 。 


表 1 A C, 序列 各 滑动 值 对 应 相关 系数 均值 


Table 1 The mean of A, Nis sliding correlation 


滑动 值 -3 d -2d -1d 0 d 1d 2d 3d 

Rai 0. 0545 0. 1072 - 0. 083 - 0.0766 -0.23 - 0.3102 -0.1512 
Rapca H5 0.0144 — 0.0923 -0. 1993 — 0.3089 —0. 3343 -0.2322 —0.204 
Raca HÀ -0. 1995 -0. 3086 — 0.4098 -0.4919 -0. 5458 -0.5576 -0. 5089 
Rot 0.3364 0.2816 0.0616 -0.1722 -0.3116 -0.37 - 0.2647 
Rape 0. 1408 0. 3873 0. 2405 - 0. 1049 - 0. 4355 - 0.4396 0.1071 
Rc, Hl -0.0515 0.217 -0.3342 0.2619 0.0247 -0.2232 0.0562 
Fac 0. 1597 0. 0528 -0.0538 -0.155 -0. 2602 -0.3653 -0.4706 
Raal 0. 1691 0. 0639 - 0.0234 - 0. 089 - 0. 1443 - 0.1953 - 0. 2493 
Ra cul —0. 0898 -0.1674 -0.2455 -0.3317 -0.4149 -0.4839 -0.5142 
Ra cat? -0.0695 -0. 1464 -0.224 -0.3079 -0.3891 -0.4583 -0.4952 
Racal! 0. 0821 0. 0984 -0.0138 -0.17 -0.3121 -0.3532 -0.2284 


X2 Efor C, 序列 各 层 信 号 取 不 同 滑动 值 时 
的 相关 系数 均值 。R_ ja,R H5S ~ Hl, 
R_j ,dL5 ~ 1A 各 为 加; Ca 原始 序列 ,高 频 5 ~ 
H1, 低频 LS ~ L1 相关 系数 均值 。 El 4 Æ foyo Ca 
序列 各 层 信号 滑动 相关 系数 统计 ,图 4( a, g) 是 
f 5 C, 原始 序列 滑动 相关 性 统计 ,原始 序列 间 表 现 
出 较 强 的 相关 性 , 且 随 滑动 值 增 大 相关 性 有 减弱 趋 
势 。 从 C, 预报 的 角度 ,我 们 希望 在 滑动 值 大 于 0 
( 滞后 相关 ) 时 县 有 更 强 的 相关 性 。 好 在 减弱 趋势 
是 缓慢 的 ,滑动 值 取 1 d 时 相关 系数 均值 为 0.65。 


—fio. 7Ca 


A(b - f) £818 f; Ci 序列 高 频 万 ~ 1 滑 
动 相关 系数 统计 ,高 频 #5, HA 部 分 较 原 始 序列 具有 
更 强 的 相关 性 , 变化 趋势 与 原始 序列 相似 。 高 频 H3 
~ Hi 部 分 则 相关 性 较 弱 , 且 变 化 不 规律 。 图 4(h ~ 
1) £888 A, , C, 序列 低频 [5 ~ L1 滑动 相关 系数 统 
计 , 低 频 部 分 相关 系数 随 滑动 值 增 大 而 减弱 ,但 整体 
表现 出 较 强 正 相 关 。 虽 然 滑动 值 小 于 0( 超前 相关 ) 
时 ,各 分 解 序列 相关 性 更 强 些 ,但 从 变化 趋势 上 来 看 
滑动 值 取 1 d 时 ,也 县 有 较 强 的 相关 性 。C, 预报 中 
可 取 滑动 值 1 d, 避 免 引 入 空间 环境 预报 误差 。 
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3 ”地磁 指数 和 阻力 系数 分 解 序列 滑动 相关 性 统计 


Fig.3 The boxplot of A, , C, sliding correlation 


X2 for Ci 序列 各 滑动 值 对 应 相关 系数 均值 


Table 2 The mean of fi, ; C, sliding correlation 


滑动 值 -3 d -24d -1d 0d 1d 2d 3d 
Ro 76 0. 639 0. 7085 0.7426 0.7224 0. 65 0.5661 0. 4805 
Ro 20, H9 0.835 0. 9464 0.9217 0. 8364 0.7526 0.6557 0.4746 
Ro 50, H4 0.705 0.7735 0. 8126 0.7866 0.7106 0.6148 0.5151 

Ryo 10,3 0.0116 - 0.0046 -0.023 — 0.0407 -0.0546 -0.0606 -0.058 
R fio. 7642 0.0018 - 0.0004 -0.0037 -0.0076 -0.0112 -0.0134 -0.0121 
R fo. 72H! 0. 0026 0. 0005 -0.0016 -0.0031 -0.0057 — 0. 009 — 0. 0106 
Ro 72415 0. 8899 0. 9958 0.8911 0.7877 0. 6868 0.5822 0.4774 

Rn 76414 0. 6288 0. 7299 0.8191 0. 8936 0. 909 0. 8634 0.774 
Ro 10,13 0. 7508 0. 8367 0. 8896 0. 9089 0. 8828 0. 8213 0. 7313 
Ro 410,12 0.728 0. 8145 0. 8491 0. 8334 0. 7751 0. 6896 0. 5899 
Ro. 56,1 0. 6809 0.7665 0.7977 0.7741 0.7 0.6054 0.5216 


上 述 分 析 表 明 阻 力 系数 C, 对 地 磁 指数 4A, 和 太 
阳 辐 射流 量 .名 ; SERE HTA IE, S A, HE fio; 
序列 的 相关 性 更 强 些 。 工 程 应 用 中 ,可 通过 实测 的 
空间 环境 指数 对 C, 序列 进行 1 ~2 d 预报 。 


3 回归 分 析 与 验证 


通过 滑动 相关 性 的 分 析 , 阻 力 系数 C, 对 地 磁 
指数 4, 和 太阳 辐射 /。; 分 别 取 滑动 值 2d 和 


1 d。 对 各 层 分 解 信 号 进行 多 元 非 线性 回归 分 
yr, 模型 为 : 


5 
C, = 2 2 
dq 7 Qo t OX; 十 QipXi + OX» t (Xp 


(6) 
采用 二 次 多 项 式 进行 回归 分 析 , 式 中 : x ,x 为 A, 
ME 分 解 的 各 层 信 号 ; Qo » Qi » (t 7 Qi3 » Aig 为 回归 
系数 ; 最 后 将 高 低频 信号 亚 加 得 到 C, (B 
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4 太阳 辐射 指数 和 阻力 系数 分 解 序列 滑动 相关 性 统计 
Fig.4 The boxplot of fi, ; ; C, sliding correlation 
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图 5 回归 分 析 结 果 统 计 与 验证 


Fig.5 The statistic of regression analysis and verification 


仍 以 解 算 的 55 d 序列 为 样本 ,选取 长 度 为 n 的 ” 4,, fi, 二 元 非 线 性 回归 模型 预报 的 C, 值 。 为 便于 
连续 序列 进行 回归 分 析 。 剩 余 序列 来 验证 回归 分 析 ” 预报 误差 比较 ,图 5( c) 依据 C, 预报 残 差 对 预报 值 
结果 ,并 统计 C, 序列 预报 残 差 。 图 5( a) 8082 的 “毛刺 ”进行 “" 谢 顶 ?平滑 处 理 。 对 比 图 5(a) 可 见 
值 分 别 取 2 d,1 d 时 4, 和 /fio; 实测 值 ,图 5(b) 给 出 “毛刺 "多 出 现在 空间 环境 上 升 或 下 降 相 。 各 回归 
不 同 回 归 模 型 预报 C 值 与 实际 解 算 值 ,曲线 标注 见 。 模型 对 空间 环境 变化 响应 不 同 ,如 f. 一 元 回归 模 
5( a) ,其 中 C, 预报 值 1 ~3 分 别 是 fi。;,4, 一 元 和 ”型 对 空间 环境 变化 处 理 较 好 ,只 含 一 个 “毛刺 ”, 其 
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余 两 模型 则 差 一 些 。 图 5(b ~ e) 均 取 n =20 d, 平 
滑 后 三 种 回归 模型 C, 序列 预报 残 差 标准 差分 别 为 : 
0. 1412 ,0. 6454,0. 2802, 可见 fio., 一 元 非 线 性 回归 
模型 预报 结果 是 最 优 。 图 5( d ~f) 给 出 fo 一 元 回 
归 模 型 时 ,分 别 取 n =30 d,20 d,10 d 时 C, 的 预报 
值 与 解 算 值 ,平滑 后 对 应 残 差 标准 差分 别 为 : 
0.1663,0.1412,0. 1580。 可 见 C, 预报 优 劣 与 回归 
分 析 序 列 长 度 相关 , 取 n =20 d 预报 结果 最 优 。 另 
外 统计 发 现 回 归 系 数 o4 ,a 0550304 波动 较 大 ， 
回归 系数 是 时 变 的 ,只 适用 于 短期 C, 预报 。 


利用 建立 的 回归 模型 选取 实例 ,检验 轨道 预报 
精度 提高 情况 。 取 空间 环境 剧烈 波动 期 样本 ,使 用 
三 种 方法 进行 预报 : O C, 取 常 数 2. 2,@ 前 一 弧 段 
估算 C, ;用 于 后 续 预 报 @) 本 文 方法 预报 C, 值 。 表 3 
列 出 三 种 方法 预报 结果 及 误差 ,可 见方 法 他 轨道 预 
报 精 度 最 优 ,位置 预报 误差 较 方法 中 、@ 的 千 米 级 ， 
精确 至 百 米 级 。 当 然 空间 环境 平稳 变化 期 的 改进 效 
果 则 没有 这 么 明显 ,限于 剧烈 扰动 期 的 数据 样本 数 
量 , 这 里 就 不 给 出 统计 结论 了 。 


表 3 不 同 C, 模式 下 的 轨道 预报 误差 


Table 3 The prediction errors of position with different C, mode 


预报 方法 C, 取 值 预报 时 长 /d 预报 X /m 预报 了 /m 预报 Z /m 位 置 误差 /m 
© 2:2 1 5941829 3075328 - 790634 2659. 573 
© 2.9304 1 5939655 3078517 - 794101 2528.256 
Q9 2.5544 1 5940775 3076875 - 792316 144. 9546 
实测 GPS 轨道 数据 5940737 3076969 - 792420 NULL 
4 ”结束语 境 变 化 段 中 表现 出 的 误差 是 不 同 的 , 故 回归 系数 对 


本 文 以 不 同 弧 长 解 算 连 续 55 d 天 宫 一 号 GPS 
轨道 数据 得 到 阻力 系数 C, 序列 ,对 同期 的 太阳 辐射 
BE fio, 和 地 磁 活 动 指数 a, A, 序列 一 并 进行 高 、 
低频 小 波 分 解 ,统计 分 析 空 间 环境 参量 4, ,fi,; 和 C, 
序列 各 分 解 信 号 滑动 相关 性 。 结 果 表 明 , 滑动 值 分 
别 取 2 d 和 1 d 时 ,卫星 阻力 系数 序列 小 波 分 解 后 
与 对 应 层次 4, ofo, 序列 具有 很 好 的 滞后 相关 性 , 验 
证 了 Ci 可 预报 性 。 地 磁 指 数 4, 和 阻力 系数 C 序列 
整体 呈现 负 相 关 , 表 明 NRLMSIS - 00 模式 在 高 度 
350 km 的 中 低 纬度 区 域 ,地 磁 活 动 上 升 相 存在 过 响 
应 ,下 降 相 则 响应 不 足 。 太 阳 辐 射流 量 fio, 和 阻力 
系数 C, 序列 整体 表现 出 较 强 正 相关 ,低频 部 分 fi, ， 
上 升 相 C, 序列 基本 保持 稳定 ,下 降 相 C, 序列 相应 
下 降 ,表明 NRLMSIS -00 模式 在 太阳 辐射 流量 上 升 
相 响 应 适度 ,下 降 相 存在 过 响应 ,高 频 部 分 呈现 很 强 
的 正 相关 ,表明 NRLMSIS - 00 模式 在 Au， 上升 相 存 
在 小 尺度 上 的 响应 不 足 , 而 下 降 相 存 在 过 响应 。 

分 析 不 同 回归 模型 和 不 同 回归 分 析 样 本 序列 长 
E n 的 情况 下 的 C, 预报 值 残 差 ,结果 表明 fio 一 元 
回归 模型 取 ”=20 d 时 , C, 预报 结果 最 优 。 另 外 统 
计 发 现 回归 系数 au,alyapyasyaa 波动 较 大 , 故 方 
法 只 适用 于 C, 短期 预报 。 此 外 磁 扰 期 C, 预报 会 出 
现 毛刺 , 需 平滑 处 理 。 因 不 同 大 气 模式 在 各 空间 环 


不 同 大 气 模式 不 具有 普遍 适用 性 。 

统计 显示 , 弧 长 3 h 解 算 阻 力 系数 序列 与 对 应 
A, 序列 呈 极 弱 相 关 性 ,可 能 是 解 算 弧 长 不 够 等 因素 
引入 其 它 误差 导致 ,具体 原因 有 待 深入 研究 。 本 文 
从 小 波 分 解 的 角度 列 出 卫星 阻力 系数 与 空间 环境 参 
量变 化 关系 ,为 C, 预报 提供 了 方法 ,解决 了 现 有 工 
程 中 普 吏 采用 前 段 定 轨 解 算 阻 力 系数 值 用 于 后 续 轨 
道 预报 时 ,在 空间 环境 剧烈 波动 期 引入 较 大 误差 的 
问题 。 分 析 方 法 对 提高 航天 器 轨道 预报 精度 具有 重 
要 参考 价值 。 
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